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Abstract. In this note we consider problems related to parabolic 
partial differential equations in geodesic metric measure spaces, 
that are equipped with a doubling measure and a Poincare inequal- 
ity. We prove a location and scale invariant Harnack inequality for 
a minimizer of a variational problem related to a doubly non-linear 
parabolic equation involving the p-Laplacian. Moreover, we prove 
the sufficiency of the Grigor'yan-Saloff-Coste theorem for general 
p > 1 in geodesic metric spaces. The approach used is strictly 
variational, and hence we are able to carry out the argument in 
the metric setting. 



I. Introduction 

The purpose of this note is to study parabolic minimizer s, which 
in the Euclidean case are related to the doubly non-linear parabolic 
equation 

3MpO-v.(iv«r^«)=o, (L1) 

where 1 < p < oo. When p = 2 we can recover the heat equation from 
(II. ip . A function u : Q x (0, T) -» R, where Q C R n is a non-empty 
open set, is called a parabolic quasiminimizer related to the equation 
( II. ip if it satisfies 

p f \u\ p ~ 2 u-^- dx dt + [ \Vu\ p dxdt 

»/supp(</)) ^ J supp(</>) 

< K I \V(u + (p)\ p dxdt, 

Jsupp(4>) 

for some K > 1 and every smooth compactly supported function in 
Q x (0, T). More precisely, in the Euclidean setting every weak solution 
to (II. ip is a parabolic minimizer, i.e., a parabolic quasiminimizer with 
K = l. 

Elliptic quasiminimizers were introduced by Giaquinta and Giusti in 
[TT| [12] . They enable the study of elliptic problems, such as the p- 
Laplace equation and p-harmonic functions, in metric measure spaces 
under the doubling property and a Poincare inequality. We refer, e.g., 
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to [3J, [5], [5], [IS]) [2D], and the references in these papers. Following 
Giaquinta-Giusti, Wieser [31] generalized the notion of quasiminimiz- 
ers to the parabolic setting in Euclidean spaces. Parabolic quasimini- 
mizers have also been studied by Zhou [321 S3], Gianazza-Vespri [TO] . 
Marchi [22], and Wang [30]. The literature for parabolic quasiminimiz- 
ers is very small compared to the elliptic case. In recent papers |18j . 
[21], parabolic quasiminimizers related to the heat equation have been 
studied in general metric measure spaces. The variational approach 
taken in these papers opens up a possibility to develop a systematic 
theory for parabolic problems in this generality. 

Our main result is a scale and location invariant Harnack inequality, 
Theorem 16.61 m geodesic metric measure spaces for a positive parabolic 
minimizer that is locally bounded away from zero and locally bounded. 
We assume the measure to be doubling and to support a (l,p)-Poincare 
inequality. We take a purely variational approach and prove the Har- 
nack inequality without making any reference to Equation (11.11) . 

In Euclidean spaces, the Harnack inequality for a positive weak solu- 
tion to Equation (11. II) . that is bounded away from zero, was proved in 
[T7] . Their proof is based on Moser's method and on an abstract lemma 
due to Bombieri and Giusti. The argument in [T7] relies on the equa- 
tion and uses, for instance, the fact that if u is a weak supersolution 
to (II. ip . then u~ l is a weak subsolution of the same equation. 

Our proof is based on the one in [17J. However, since we deal with 
parabolic minimizers and upper gradients in the metric setting, changes 
in the argument are required. To give an example, in the strictly vari- 
ational setting it is not true that if u is a parabolic super minimizer, 
then u~ l is a parabolic subminimizer. Instead we establish the required 
estimates separately for both super- and subminimizers. 

Grigor'yan [13J and Saloff-Coste [25J observed independently that 
the doubling property and a Poincare inequality for the measure are 
sufficient and necessary conditions for a scale and location invariant 
parabolic Harnack inequality for solutions to the heat equation (p = 2) 
on Riemannian manifolds. Later Sturm [29] generalized this result to 
the setting of Dirichlet spaces. 

One motivation for the present note is to show the sufficiency for 
general 1 < p < oo in geodesic metric measure spaces without invoking 
Dirichlet spaces or the Cheeger derivative structure for which we refer 
to [H]. We also refer to a recent paper [2] and to [T] on parabolic 
Harnack inequalities on metric measure spaces with a local regular 
Dirichlet form. It would be very interesting to know whether also the 
necessity holds in this general setting. 

Very recently a similar question has been studied for degenerate par- 
abolic quasilinear partial differential equations in the subelliptic case by 
Caponga-Citti-Rea [8] . Their motivating example is a class of subellip- 
tic operators associated to a family of Hormander vector fields and their 
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Carnot-Caratheodory distance. The setup in the present paper cover 
also Carnot groups and more general Carnot-Caratheodory spaces. 



2. Prelimininaries 

In this section we briefly recall the basic definitions and collect some 
results we will need in the sequel. For a more detailed treatment we 
refer, for instance, to a monograph by A. and J. Bjorn [I] and to 
Heinonen [T3], and the references therein. 



2.1. Metric measure spaces. Standing assumptions in this paper 
are as follows. By the triplet (X, d, /i) we denote a complete geodesic 
metric space X, where d is the metric and fi a Borel measure on X. 
The measure fi is supposed to be doubling, i.e., there exists a constant 
Cji > 1 such that 

< fi{B(x, 2r)) < C^{B(x, r)) < oo (2.1) 

for every r > and Here B(x, r) := {y G X : d(y, x) < r}. We 

denote \B = B(x, Ar) for each A > 0. We want to mention in passing 
that to require the measure of every ball in X to be positive and finite is 
anything but restrictive; it does not rule out any interesting measures. 
Equivalent ly, for any x G X, we have 

rtB(«,fl)) /fly (2.2) 



n(B(x,r)) 

for all < r < R with g M := log 2 C M , and C > is a constant which 
depends only on C^. The choice = log 2 is not necessarily optimal; 
the exponent q^ serves as a counterpart in metric measure space to the 
dimension of a Euclidean space. In addition to the doubling property, 
we assume that X supports a weak (l,p)-Poincare inequality, see below. 
Moreover, the product measure in the space X x (0,T), T > 0, is 
denoted by v = fi ® C 1 , where C 1 is the one dimensional Lebesgue 
measure. 

It is worth noting that our abstract setting causes some, perhaps 
unexpected, difficulties. For instance, in not too pathological metric 
spaces it may happen that B(x\,r\) C B(x2,r 2 ) but B(x2,2r 2 ) C 
B(x 1 ,2r l ). 

We follow Heinonen and Koskela [15] in introducing upper gradients 
as follows. A Borel function g : X — > [0, oo] is said to be an upper 
gradient for an extended real-valued function u on X if for all paths 
7 : [0, lry\ — > X, we have 

|u(7(0))-u(7(* 7 ))l < / 9ds. (2.3) 

If (12.31) holds for p-almost every path in the sense of Definition 2.1 
in Shanmugalingam [27], we say that g is a p-weak upper gradient of u. 
From the definition, it follows immediately that if g is a p-weak upper 
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gradient for u, then g is a p-weak upper gradient also for u — k, and 
\k\g for ku, for any fcGl. 

The p-weak upper gradients were introduced in Koskela-MacManus 
[2Tj . They also showed that if g G L P (X) is a p-weak upper gradient 
of u, then one can find a sequence of upper gradients of u such 

that gj —7- (7 in L P (X). If w has an upper gradient in L P (X), then it has 
a minimal p -weak upper gradient g u G in the sense that for every 

p-weak upper gradient g G L P (X) of u, g u < g a.e., see Corollary 3.7 
in Shanmugalingam [28] . 

Let Q be an open subset of X and 1 < p < 00. Following the 
definition of Shanmugalingam [27] . we define for u G L p (f2), 

IHIiv 1 ^) = + llfl''"llxp(n)' 

where the infimum is taken over all upper gradients of u. The Newto- 
nian space iV 1,p (f2) (c L p (fl)) is the quotient space 

jV^O) = \u G Z7(ft) : ||u||jvi,p (n) < 00} / ~, 

where u ~ u if and only if ||« — ulljvLj^n) — 0. The space N 1,P (Q) is a 
Banach space and a lattice, see Shanmugalingam [27]. If w, f G A^ 1,P (J1) 
and u = v /i-a.e., then u ~ t> . However, if w G A^ 1,p (u;), then u ~ u if 
and only if -u = f outside a set of zero Sobolev p-capacity [27] . 

A function w belongs to the /oca/ Newtonian space N^(Q) if w G 
7V 1,P (V) for all bounded open sets V with K C Q, the latter space 
being defined by considering V as a metric space with the metric d and 
the measure \i restricted to it. 

Newtonian spaces share many properties of the classical Sobolev 
spaces. For example, if u, v G N^(Q), then g u = g v a.e. in {x G 
tt : u(x) = v(x)}, in particular g min{tiiC} = g u X{u^c} for c G R. 

Remark 2.4. Note that as a consequence of the definition, the functions 
in N 1,P (Q) are absolutely continuous on p-almost every path. This 
means that u o 7 is absolutely continuous on [0, length^)] for p-almost 
every rectifiable arc-length parametrized path 7 in Q. This in turn 
implies that for each of these paths we have \(u o 7)'(s)| < 0(7(5)) for 
almost every s G [0, length^)]. We refer to Theorem 1.56 and Lemma 
2.14 in [1]. 

We shall also need a Newtonian space with zero boundary values. 
For a measurable set E C X, let 

N^iE) = {f\ E : / G N X >P(X) and / = on X \ E}. 

This space equipped with the norm inherited from N 1,P (X) is a Banach 
space. 

We say that X supports a weak (l,p)-Poincare inequality if there 
exist constants C p > and A > 1 such that for all balls B(xq,t) C X, 
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all integrable functions u on X and all upper gradients g of u, 



4- \u — Ub\ dfx < C p r I 4 g p d[i\ , (2.5) 

J B(x ,r) \J B(x ,Ar) J 



where 



u B ■= f u ^'- = ~Tnt vT / ud f J '- 

'B(xo,r) V>{-E>{X0, r)j JB{x ,r) 



If the metric measure space X has not "enough" rectifiable paths, it 
may happen that the continuous embedding iV 1,p — > L p , given by the 
identity map, is onto. If X has no nonconstant rectifiable paths, then 
g u = is the minimal p-weak upper gradient of every function, and 
N 1,P (X) = L P (X) isometrically. The fact that the Newtonian space is 
not simply L P (X) is implied, for instance, by assuming that X supports 
a weak (l,p)— Poincare inequality. 



2.2. Parabolic setting. Our set-up is the following. Let f2 C X be 

an open set, and < T < oo. We write Qt '■= ^(o,t) := ^ x (0, T") 
for a space-time cylinder, and z = (x,t) is a point in Q?. We denote 
by L p (0,T; N l,p (Q)) the parabolic space of functions u : Qt K such 
that for a.e. t G (0,T), x >-)■ u(x,t) belongs to A^ 1,p (f2) and 

T 

and similarly for Lf oc (0,T; N 1( £(Q)). Here we have defined 

</«(ar,t) := 

at //-almost every (x,t) G x (0,T). 

The following calculus rules will be used throughout the text. As- 
sume u, v G Ll c {0,T;N^{n)). Then for almost every t and //-almost 
every x 

gu+v < 9u + gv, 

9uv < \u\g v + |u|sr u . 

In particular if c is a constant, then g cu = \c\g u . For proof at each time 
level, see [I]. This proof guarantees that g u+v and g uv are defined at 
almost every t and //-almost every x. The definition of the parabolic 
minimal p-weak upper gradient then implies the result. Note that the 
above does not claim that uv is in the parabolic Newtonian space, even 
if u and v are. 

In the Euclidean case it can be shown that stating that a function 
ti:ftx(0,r)4l,!ie £L(°> T ; W ioc( n )) is a weak solution to the 
doubly nonlinear parabolic equation (11. ip . is equivalent to stating that 
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u is fulfills the variational problem 

p [ \u\ p ~ 2 u^- dx dt + [ \Wu\ p dxdt 

J supp(</>) J supp(</>) 

< / \Vu + V<p\ p dxdt, 

J supp(0) 

for every <fi G C^°(Q x (0, T)). Since partial derivatives cannot be 
defined in a general metric space, there is little sense in trying to define 
the weak formulation of Equation (II. ip in the metric setting. The 
variational approach on the other hand only considers integrals with 
absolute values of partial derivatives and an inequality - as opposed 
to demanding a strict equation with gradients. This opens up the 
possibility to extend the definition of a parabolic minimizer related to 
the doubly nonlinear equation to metric measure spaces in the following 
way: 

Definition 2.6. We say that a function u G £f oc (0, T; is a 

parabolic minimizer if we have 

pf \u\ p - 2 u^du+ [ g p u dv<\ g p u+<p du, (2.7) 

J supp(0) J supp(<j!>) J supp(<?!)) 

for all G Lip (fi r ) = {/ G Lip(fi r ) : supp(/) C Vl T }. If (J23D holds 
for all nonnegative G Lip (f2T) a function it G £f oc (0, T; N^^l)) is a 
parabolic superminimizer; and a parabolic subminimizer if ( 12. 7p holds 
for all nonpositive G Lip (fir)- 

Observe that here parabolic minimizers are scale invariant but not 
translation invariant. 

Let < a < 1, let parameters r and T be positive, and to G R. A 
space-time cylinder in X x R is denoted by 



Q ar {x, t) = B(x, or) x (t - T(ar) p , t + T(ar) p ). 
It will also be of use to define positive and negative space-time cylinders 



as 



1 — a\ p fl + a 

aQ + (x, t) = B(x, ar)x ( t + T [ — — J r p ,t + T 



l + a\ p /l-a XP 
a Q-(x,t) = B{x,ar) x [t-T[ — — - r p ,t-T' 



■p 



Using these, we write 

Qar = Qar(x ,t ), OlQ + = CiQ + (x ,t ) , aQ~ = aQ~ (x , t Q ) . 

Above r is chosen according to (xq, to) and T in such a way that Q r C 
Qt- Our goal in this note is to prove the following Harnack inequality 
using Moser's argument and energy methods: 
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Suppose 1 < p < oo and assume that the measure \i in a geodesic 
metric space X is doubling with doubling constant C^, and supports a 
weak (l,p)-Poincare inequality with constants C v and A. Then a par- 
abolic Harnack inequality is valid as follows: Let u > be a parabolic 
minimizer in Q r C Qt, locally bounded and locally bounded away from 
zero. Let < 5 < 1. We have 

ess sup u < C ess inf u, (2.8) 

6Q- «2+ 

where C = C(C M , C p , A, T). 

Note that the constant in the Harnack estimate does not depend on 
r and so is scale invariant, as long as r is such that Q r C Qt- The 
parameter T controls only the relative proportions of the spatial and 
time faces of Q r . 

2.3. Sobolev— Poincare inequalities. We shall need the Sobolev in- 
equality for functions with zero boundary values; if / G Nq P (B(x , R)), 
then there exists a constant C > only depending on p, C M , and the 
constants C p and A in the Poincare inequality, such that 

( / 7 <Ctt( / g p f dfi) '\ (2.9) 

\Jb(x ,R) J \JB(xo,R) J 

where 

( oo, if p > g M . 

For this result we refer to [20]. The following weighted version of the 
Poincare inequality will also be needed. 

Lemma 2.10. Let f e N^ p (B(x , R)), and 

d(x, xq) 



4>{x) = l - 



R 



B(xo,r) J B(xo,r) 

for all < r < R, where 



where 9 > 0. Then there exists a positive constant C = C(Cn,C p ,p, 
such that 

I |/-^| p 0^<Cr p / g^dfx (2.11) 

J B(xn,r) J B(xn.r) 

* Ibm^^' 

Sketch of proof. The main idea in the proof, for which we refer to Saloff- 
Coste [26J Theorem 5.3.4], is to connect two points in the ball B(xo,r) 
with a certain finite chain of balls. For this chain we need to assume 
that our space X is geodesic. □ 

7 



Remark 2.12. Lemma ^. 101 will be used later in the proof of Lemma I5T21 
We stress, however, that apart from Lemma 12.101 all other estimates 
prior to Lemma 15.21 are valid without X being geodesic. 

2.4. Bombieri's and Giusti's abstract lemma. A delicate step in 
the proof based on Moser's work is to use a parabolic version of the 
John-Nirenberg inequality, i.e., exponential integrability of BMO func- 
tions. To avoid the use of the parabolic BMO class, the parabolic 
John-Nirenberg theorem is replaced with an abstract lemma due to 
Bombieri and Giusti [7]. Consult [26J or [17J for the proof. 

Let v be a Borel measure and consider a collection of bounded mea- 
surable sets U a , < a < 1, with U a > C U a if a' < a. 

Lemma 2.13. Fix < 5 < 1, let #,7, and A be positive constants, 
and < q < 00. Moreover, if q < 00, we assume that 

v{U x ) < Au(U s ) 

holds. Let f be a positive measurable function on U\ such that for every 
< s < min(l,g/2) we have 

for every a, a' such that 0<5<a'<a<l. Assume further that f 
satisfies 

log/>A})<^ 

for all A > 0. Then 




where C = C(q,5,9,j,A). 

3. Reverse Holder inequalities for parabolic 
superminimizers 

In this section, we prove an energy estimate for parabolic supermini- 
mizers. After this, using the energy estimate we prove a reverse Holder 
inequality for negative powers of parabolic superminimizers. 

Establishing energy estimates for parabolic superminimizers is based 
on substituting a suitably chosen test function into inequality (12. 7p . 
and then performing partial integration to extract the desired inequal- 
ity from it. While doing this, we take the time derivative of w p_1 , 
even though u is not assumed to have sufficient time regularity for 
this. Therefore, the reader should consider the time derivation of u as 
being formal. Justifications for the formal treatment will be given in 
Remark 13.61 
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Lemma 3.1. Let u > be a parabolic superminimizer, locally bounded 
away from zero, and < e ^ p — 1 . Then 



ess sup 
o<t<T Jn 



u 



p-l-ein 



p dfi + 



' supp(4>) 



sftur 1 -*^ dv 

+ C 2 [ u p - l - e \{cf?) t \dv, 



' supp(4>) 

for every <f> G Lip (f2 T ), < < 1, where 
e\p — 1 — e 



supp(<f>) 



p\p 

e 



1 + 



Pip 



Co 



1 + 



Kg - 1) 

e\p- 1 - 



Proof. Assume e > 0, e ^ p — 1. Let be a function < < 1, 
4> G Lip (f2;r). Since has compact support, we can choose < t\ < 
t 2 < T such that <f)(x,t) = yU-almost everywhere when t G" (ti,t 2 ). 
Since by assumption u is locally bounded away from zero, we may 
assume a positive constant a > such that after denoting v = au we 
have 1 — ecf^v' 6 ^ 1 > //-almost everywhere in the support of 0. It then 
follows that //-almost everywhere in the support of 0, we have 



9v+<j>Pv- 



<{X-e4Pv- e - x )g v +pv- e 4T%- 



(3.2) 



That (13. 2p does indeed hold can be seen the following way: Let t be 
such that v(-,t) G N 1,P (Q). Consider any arc-length parametrization 
7 of a rectifiable path on which v(-,t) is absolutely continuous. Since 
</>(•, t) is Lipschitz-continuous, it is absolutely continuous on 7 . Define 
h: [0,length( 7 )] -> [0,oo), 

h(s) = u( 7 (s), t) + 0(7(5), t) p i;(7(s), t)- £ . 
Then ft, is absolutely continuous, and so we have 



h'(s) 



+ P <l>( 1 {s),ty- 1 v( 1 {s),ty 



ds 

ds 



for almost every s G [0, length( 7 )] with respect to the Lebesgue mea- 
sure. We know that 



dv(j(s),t) 



ds 



< g v (j{s),t), 



ds 



< 94>(l{?), t), 



for almost every s G [0, length( 7 )]. Hence 



d(v + 4> p v- £ )(j(s),t) 



ds 



= \h'(s)\ 

< (1 - e0( 7 (s), tyv( 7 (s), ty^gMs), t) 
+ p<j ) ( 1 (s),ty- 1 v(l(s),tr e g^ 1 (s),t), 



for almost every s G [0, length^)]. The fact that this holds for p-almost 
every rectifiable path 7 now implies (13. 2p . Using the convexity of the 
mapping it K f we have 



a p 

•3v-\-<f)Pv~ 



< (1 _ ^^-1)^ + ^1-^-1^. 
Assume < T\ < r 2 < T. For a small enough h > 0, denote 



Xfn,«](*) 



' (t - Ti)/h, Ti<t<t<T x + h 

1, r 1 + h<t<r 2 -h 

(r 2 - £)/7i, r 2 -h<t <r 2 

0, otherwise. 



(3.3) 



Integrating by parts, we find 



(P-1) 

p — 1 — £ 



Tl ./ft 



After taking the limit /i — >■ in the expression above, and using 
Lebesgue's differentiation theorem, we have 



lim 



r-2 



,p-l- £ ^p 



(i/i 



Qx{t} 



+ 



t=Ti 



p — 1 — £ 



v v - x - e {^P) t dndt 



As m is a positive parabolic superminimizer related to the doubly non- 
linear equation, also v is a parabolic superminimizer. Moreover, by 
remark 13.61 below, <fr p v ~ e x\ Tl T2 i is a nonnegative admissible test func- 
tion. Hence by the definition of a parabolic superminimizer and (13. 3 p 
we have 



Pip ~ 1) 
p — 1 — e 



< lim 



< -£ 



dfi 



JUx{t} 



T2 



+ 



t=Tl 



T2 



Tl JSi 



supp(</-Pxf TliT2] ) 



9 P v+4>Vv -e x h dv 



£ - 1 gPdu + p p £ 1 - p 



SUpp(<^X[ri,T 2 ]) 



?- l - e qldv. 



supp(^X[ n ,T 2 ]) 
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On one hand, setting t\ = t\, and r 2 = t 2} we obtain 



Jti Jn e\p-l- e\ J tl J n 




\(0 p ) t \dfidt 



(3.4) 



9% d(i dt. 



On the other hand, if e < p — 1, set T\ = t and t 2 = t 2 . If e > p — 1, 
set ri = t\ and r 2 = t. We obtain 



This holds for almost every t G (t 1 ,t 2 ). Dividing (13. 5p by p{p— l)/e\p — 
1 — e|, and adding the resulting expression to (13. 4p yields the desired 
estimate for v, since the constants in the inequality do not depend 
on i 6 (^1,^2) and 0, g$ vanish outside the support of 0. The proof 
is completed by dividing the resulting expression sidewise with the 
constant a p ~ 1 ~ £ . □ 

Remark 3.6. We now give justifications for the formal treatment above. 
By a change of variable, it is straightforward to see that for a nonneg- 
ative parabolic super- or subminimizer v and for an admissible test 
function ip, for any small enough s, we have 



We multiply this inequality sidewise with a standard mollifier with re- 
spect to the time variable s, and then integrate both sides with respect 
to s. After using Fubini's theorem on the left side, this yields 




(3.5) 





(3.7) 



Here we have used the notation 
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where 8 is a standard mollifier and a > is assumed to be small enough 
so that everything stays in the time cylinder. To be precise, in the proof 
of Lemma 13.11 we then choose the test function 



with GLip (O( tlj t 2 )). The test function ip now has compact support 
and belongs to the space L p (0, T; N l,p (Q)). By Lemma 2.7 in [23], easily 
adaptable for minimizers related to the doubly nonlinear equation, ip 
can be plugged into inequality (13. 7p . Similarly to the formal proof 
above, partial integration is then performed to write the expression in 
a form where (v p ~ l ) a is not differentiated with respect to time. Once 
this is done we can take the limits o — > and h — > 0, which leads 
us back to the inequality above (I3.4p . For details on justifying the 
convergence of the upper gradient terms in (13. 7p as a — > 0, we refer the 
reader to [2"3] . 

Next we prove a reverse Holder type inequality for negative powers 
of parabolic superminimizers. The first step of the proof consists of 
combining Sobolev's inequality with the energy esimate of Lemma 13.11 
Then, because the energy estimate is homogeneous in powers, the ob- 
tained inequalities can be combined as in Moser's iteration to complete 
the proof. 

Lemma 3.9. Let u > be a parabolic superminimizer in Q r C Qt, 
locally bounded away from zero and let < 5 < 1 . Then there exist 
constants C = G[G^ C p , A,p, 5, T) and 9 = 8(0^, p) such that 



for every 0<5<a / <Q;<l and for all < q < p. 

Proof. Let us fix a', a such that 0<5<a'<a<l, and divide the 
interval (a', a) as follows: cto = a, ctoo = a', and 



Qj = Q aj r = Bj x Tj = B(x , af) x (t - T(a j r) p , t + T(a j r) p ), 



and choose the sequence of test-functions {0j}°^ o so that supp(^) C 
Qj, < <pj < 1 on Qj, and 4>j = 1 in Qj+i- Moreover, let each <pj be 
such that 



V> = ^((Oa)^X[U]> 



(3.8) 




ctj = a — (a — a')(l — 7 J ), 
where 7 = 2 — p/n = I + (k — p) /k> 1. We set 





\p - 1 - £ 



) 



P 



P 



U 
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Using Holder's inequality brings us to the following estimate 
u^ 1 -^ dfldt 



3 + 1 " B j + 1 



{k- P )/k / \ p/K 

<4 [4 u^-^tfdfi) \ J- {u^~ l ~ £)/p (t) 3 ) K dA dt 



T j+1 B j+1 j \j a J+1 

(\ (k—p)/k 
ess sup / u^'^dp 
Tj J Bj J 

■j (j (u {p - 1 ~ £)/p (l) j ) K df)j dt. 



Observe that \Tj\ = 2T(ajr) p and <x, + i > min{5, (1 + 7) _1 }<x,-. Thus 
the multiplicative constant on the right hand side is bounded by a 
constant independent of j, r, T, a', and a. We estimate the last term in 
the preceding inequality by Sobolev's inequality, see (12. 9p . We find 



u^-^r dp ) < Cr p £ < ( ,_ 1 _ e)/p ^ dp 



< Cr p ^ [vP- X -% + ( fe * £] ) u-^g p J p ) dp, 



where C = C{C^,C P , K,p). Since e > 0, e ^ p — 1, we may use 
Lemma [3. II to obtain 

J- J- u^- 1 -^ dpdt<cl ess sup -I u p ~ l - e (t) p dp j 



(k-p)/k 

< c 1 1 1 ' ' x 



C 

T5p 



T 3 J B 3 



u^fc^. + C^^dpdt 

+ b - 1 - elV- 1 -*^. + C 2 |(#) t |)) dp eft 



where C = C(e,C^,C p ,A,p,5,T) is uniformly bounded for every e, 
except in the neighborhood of e = 0. For each j = 0,1, . . . we can now 
use the above with Ej > 2p — 1 chosen in such a way that p — 1 — Sj = 

13 



P7"'; to write 



u- p ^ +1 dv 



■j+l 



jp \ -l/pyi / r \ -l/p7 J 



(3.10) 



where C = C(C^, C p , A,p, 5, T). By iterating this, since 7 > 1, we find 
that 



ess inf w > 

v^oo -1 v^ 00 -(! + 7)j / I \ ^J^ ~ 



[a — a') 



C \ ( f _ \ - 1/ p 
u p dv 



(a — a 



0) (£ 



where the constant C = C(C^, C p , A,p, 6, T) is positive and finite. The 
proof is now completed for any < q < p by using a result from real 
analysis, see Theorem 3.38 in [IB]. □ 



We also prove a reverse Holder inequality for positive powers of par- 
abolic superminimizers. 

Lemma 3.11. Let u > be a parabolic superminimizer in Q r C Qt 
which is locally bounded away from zero, and < 5 < 1 . Then there 
exist constants < C = 0(0^, C p , A,p, q, S, T) and 6 = 9(C^,p) such 
that 

for all < 5 < a' < a < 1 and for all < s < q < (p - 1)(2 - p/k) 
and k is as in (12.31) . 



Proof. Assume < s < q < (p — 1)(2 — p/k), where k is as in the 
Sobolev-Poincare inequality. Then there exists a k such that s^ k ~ l < 
q < s^ k . Let po be such that < po < s and q = 7 fc po- Now for each 
j — 0, k — 1, there exists a < Sj < p — 1 such that p — l — ej = po7 J . 
By the first part of the proof of the previous lemma, we have 



u p0 ^ +1 dv 



3+1 



i/poy 



2+1 



,,p:i 



i/po7 j / r \ 1 /poy 
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where C = C(C^ C p , A,p, q, 8, T). Iterating this estimate for j = 
0, . . . , k — 1 yields 

< 7 rr) f u po dv\ , (3.12) 

where C = C(C M , C p , A,p, q, 8, T) blows up as q tends to (p~l)(2—p/ k) 
and 

7 . = j^ ( i- 7 -)<^. 

7 — 1 7 — 1 

Using Holder's inequality on the right hand side of (I3.12p . setting 9 = 
VI l(l — 1) an d using the fact that s/j < po < s completes the proof. 

□ 



4. Reverse Holder inequalities for parabolic 
subminimizers 

In this section we prove estimates analogous to those in Section [31 
but this time for parabolic subminimizers. This is done essentially 
identically to what was done for superminimizers, but with a slight 
change in the test function we use. Then we utilize the obtained energy 
estimate to prove a reverse Holder inequality for positive powers of 
parabolic subminimizers. 

Lemma 4.1. Let u > be a locally bounded parabolic subminimizer 
and let e > 1 . Then 

ess sup / u p ~ 1+£ (f) p dp + / u £ ~ 1 g u (j) p dv 

0<t<T Jq ' Jsupp(<l>) 

<C X \ u p - 1+£ g p dv + C 2 [ u p ~ 1+£ \(r)t\ dv, 

J supp(<f>) J supp((f>) 

for every <ft 6 Lip (QT), < <fi < 1, where 

Proof. Let < < 1, G Lip (Q T ), for some < ti < t 2 < T. Let 
e > 0. Since by assumption u is locally bounded, we can take a constant 
a > such that after denoting v = au, we have 1 — ecjfv^ 1 > almost 
everywhere in the support of 0. Since u is a subminimizer, also v is a 
subminimizer and we can plug — 0(x, t) p v(x, t) £ test function 

into inequality (12. 7p . The rest of the proof is now completely analogous 
to the proof of Lemma 13.11 □ 

Now we prove a reverse Holder type inequality for positive powers 
of parabolic subminimizers. Again, the proof consists of combining 
the energy estimate of Lemma T4. II with Moser's iteration to obtain the 
inequality. 
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Lemma 4.2. Let u > be a parabolic subminimizer in Q r C 
which is locally bounded and let < 5 < 1 . Then there exist constants 
C = C(C^,C p ,A,p,5,T) and 9 = 0{C^,p) such that 

( C \ 1/q f f \ 1/q 

ess sup u < ( -, r~r I I + u q dv\ 

holds for every 0<5<a'<a<l and for all < q < p. 

Proof. The steps of the proof are analogous to the proof of Lemma 13.91 
The difference is that here we use Lemma 14.11 and the observation that 
for each 7 J , j — 0, 1, . . . there exists a Sj > 1 such that p — l + ej = p^ . 

□ 



5. Measure estimates for parabolic superminimizers 

The following logarithmic energy estimate will also be important 
to our argument. Regarding the time derivation of m p_1 , the proof 
presented below is again formal. Justifications for this can be given as 
in Remark 13.61 we use the test function as in (13. 8p . but with e = p— 1. 

Lemma 5.1. Let u > be a parabolic superminimizer, locally bounded 
away from zero. Then 
" T2 

gfogu^ d ^dt - p 



T2 

log u<^ dfi 

_Jnx{t} 



t=Tl 



< f 



( P - 1)P 



/•T2 r 

g^d/j,dt + p / |logw||(0 p ) t | dfidt 



holds for every <f> G Lip (QT), such that < < 1 and almost every 

o < n < r 2 < T. 

Proof. Let (ft G Lip (f2<r) be such that < <fi < 1. As in the preceding 
lemma, since both the definition of a parabolic superminimizer and the 
final weak Harnack inequality are scalable properties, we may assume 
that u has been scaled in such a way that 1 — (p — l)<f) p u~ p > almost 
everywhere in the support of <f), and by using the convexity of the 
mapping t i-> t p , we find 

< (l - (p - i)<fu- p )g p u + p p (p - i) 1 ^- 

Let X[ n T2 ] be defined as in Lemma 13.11 Integrating by parts, we obtain 



T2 



u 



O-'^U^X^tdfldt 



(p - 1) / lJ\oguW(X[ T1 ,r 2] )t + (log«)(^) t X[n.«]] dlidt. 
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Taking the limit h — > 0, we obtain by Lebesgue's theorem of differen- 
tiation 



lim r j u^u^ X \ TiM ) t d^dt 



= (P-1) " 



(log u) 4> p d/i 



Jix{(} 



+ 



t=Tl 



ri Jn 



(logw)(0 p )t <i/i<it 



As u is a parabolic superminimizer and <fi p u l p x\ Tl T2 \ is a non-negative 
admissible test-function, we obtain 



p(p - 1) - 



log d/x 



f2x{t} 



+ 



T2 



(log«)(0 P M/idt 



< lim 



supp(</>Pxf TliT2 j) 

T2 



^supp(^ X f T1 , T2] ) [T1 ' T2] 



<-(p-l) / / ^ u d^dt + ^(p-l) 
Rearranging terms completes the proof. 



i-p 



n </n 



| d/i (it. 



□ 



Next, using the logarithmic energy estimate, we establish monotonic- 
ity in time of the weighted integral of logw. This in turn enables us to 
estimate the measure of the level sets of log-u around a time level to- 

Lemma 5.2. Let u > be a parabolic superminimizer in Q r C Qt and 
assume u is locally bounded away from zero. Let < a < 1. Define 

where < a < 1 and (x, t) G Q r . Let 

/? = i-/ logM(x,t o )0 p (a;)rf/i, 

N — f <P P (x)dfi. 

JB(x ,r) 

Then there exist positive constants C = C(C^,C p ,p,a) and C = 
C'(C^,p,a) such that 

MQar) 



V 



({(x,t)eQ ar : t<t ,\ogu(x,t)>\ + p + C'})<c'^^ 



and 



v({{x,t)eQ ar : t>t ,\ogu(x,t)<-\ + f3-C'}) <C^^ 
for every A > 0. 
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Proof. From the definition of it readily follows that < <f> < 1, 
g<t> < (ar) _1 , and for every t E [to — T(ar) p ,t + T(ar) p ], 

We write 

1 f 

v (x, t) = logw(x, t) — f3 and V(t) = — / v[x, t)4> p (x) dfi, 

N Jb(x ,t) 

and find that V(to) = 0. Let < £(t) < 1 be a smooth function 
such that supp(0 C (t - Tr p ,t + Tr p ), and f(t) = 1 for all t E 
[t -T(ar) p ,t + T(ar) p ]. 

Write ip(x, t) = <p(x)C,(t). Since u is a positive superminimizer bounded 
away from zero, we can use Lemma 15.11 with ip as a test function. We 
obtain for t — T(ar) p < ti < t 2 < t + T(ar) p , since on this interval 

m = i, 

VI gl<fdiidt-p\NV(t)\ t2 <-^—-[ t2 [ g p d^dt 

Jtr JB(x ,r) L J *=*i [p-lFJt! JB(x ,r) 



'h JB(xo,r) 

Cjf , KB(x ,r)) 

-J^iy (2_l) (ar)P ' 

where C = C(p). On the other hand, from the weighted Poincare 
inequality ( 12. lip , we have 



2 



\v -V(t)\ p d^dt 

ti J B(xQ,ar) 

< I' [ \v-V{t)\ p (tfd^dt<Cr p [ 2 [ g p (f) p dfidt } 

where C = C(C M , C p ,p, a). By combining these we find 

1_ " r /_ ' \v -Vit^dfidt + Vih) -V(t 2 ) 



CNr p 



tl J B(xo,ar) 



< C(t 2 -ti)//(^(x ,r)) < ( 2 \ p C{t 2 -t x )ii{B{x Q ,r)) 



N(ar) p \1 — a J jj,(B(x , ar))(arY 

<c , (t2-h) 

(ar) p ' 

where C = C'{G^ C p ,p, a). We denote 

/ \ / \ C (t — to) 

w{x, t) = v(x, t) + — — — 

{ar) p 

and 

C'(t - 1 ) 
(ar) p ' 
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W(t) = V(t) + 



and restate the preceding inequality as follows 
(1 - a) p r t2 



CNrP 



I 2 I \w- W(t) \ p dfidt + W(h) - W(t 2 ) < 0. 

Jtj J B(x ,ar) 



\w -W(t)\ p dfi-W'(t) < (5.4) 



This implies that W{t x ) < W(t 2 ) whenever t - T(ar) p < h < t 2 < 
to + T(ar) p , i.e. the function W is increasing, thus differentiable for 
almost every t G (to — T(ar) p , t Q + T(ar) p ). As a consequence we obtain 

CNr p J Bixo , ar) 

for almost every to — T(ar) p < t < t + T(ar) p . Let us denote 
Exit) = {x G B(x ,ar) : w(x,t) > A}, 

E\ = {( X it) ^ Qar '■ t < t , W(x,t) > A}. 

For every t - T(ar) p < t < t and A > 0, since W(t) < W(t ) = 0, we 
have 



(A - W{t)Y^E-{t)) < \w- W(t)\ p d^. 

' B(xQ,ar) 



Hence we have 



%(£-(t))- - <0, 



CNr p ^ v A w/ (A - W(t))p 

for almost every t — T(ar) p < t < t . This yields, after integrating 
over the interval (t — T(ar) p , t ), 

^ < TT^P - W(t))-^]l tM < ° 



NrP ~ (1 - v y; k=t -T( a r)p - (i _ a^AP- 1 

where C = C(C^ C p ,p, a). Together with (15. 3p . this implies 



^{(M)eQ ar : t<t , hgu(x J t)>X + /3 + C'}') <C^pp, 

where C = C(C M , C p ,p, a). Denote then 

E^(t) = {x G B(a;o, ar) : w(x,t) < -A}, 
£+ = {( X) t) eg ar : t > t , w (x,t) < -A}. 

Similarly to the case of E^~, using the monotonicity of W(t), we obtain 



(\ + W(t)) p ^E+(t))< / \w-W(t)\ p dfi, 

J B(xo,ar) 

for every to < t < to + T(ar) p . This together with (15.41) leads to 
(l-an(E+(t)) W'[t) 



CNr p (A + W(t))P 
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< 



for almost every i < t < to + T(ar) p . Integration over the interval 
{to, t + T(ar) p ) gives now 

iVr p " (1 - a)» K + WJ Jt=t0 " (1 - ay\v-i ' 
and thus after using (15. 3p we may conclude 

v({{x, t)EQ ar : t> to, \ogu<-\ + /3- C'}) < C U -^-. 
Again C = C(C^, C p ,p, a). □ 



6. Harnack's inequality for parabolic minimizers 

Having established a logarithmic measure estimate for superminimiz- 
ers around a time level to, we have the prerequisites to use Lemma f2. 131 
This way for parabolic super minimizers we can glue the reverse Holder 
inequality for negative powers together with the reverse Holder inequal- 
ity for positive powers. We obtain a weak form of the Harnack inequal- 
ity for parabolic superminimizers locally bounded away from zero. This 
result is in some sense finer than the final Harnack inequality since we 
only assume the superminimizing property, and hence it is of interest 
in itself. Observe in the following how, from applying Lemma 12.131 
separately on both sides of the time level to, a waiting time inevitably 
appears between the negative and positive time segments. 

Lemma 6.1. Let u > be a parabolic superminimizer in Q r C Qt 
which is bounded away from zero. Then 



1/9 

< C ess inf u, 

<6Q- J ^ 



where < 5 < 1 and < q < (p — 1)(2 — p/k). Here 
C = C(C f ,,C p ,A,p,q,6,T). 

Proof. Assume < 5 < 1. Let /3 and C be as in Lemma 15.21 By 
Lemma 1331 there exists a positive constant C = C{C P ,C P , A,p,S,T), 
such that for every < s < p and 0<<5<o; / <a;<l,we have 



(ess sup u 1 e 13 °') l = ess inf ue /3+c " 

a'Q+ a'Q + 
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-x/s (6.2) 



By Lemma [5.21 applied to { Q3+s r : t > t }, we have 

5 



*({0M) e = log(«" 1 c /, "°') > ^ 

< i/({(x,t) G Qs±s r : t>t , logfa-V" ') > A}) (6.3) 



Ap-i - Ap- 1 ' 
for every A > 0. In the last step of the above inequality we used the 
doubling property of fi, and so C = C(C^,C p ,p,5). From f 16 . 2 j) and 
f )6.3p . we now see that the conditions of Lemma l'2.1'd\ with ^-Q + in 
place of Ui, are met. Hence 

esssuptrV- " < C, (6.4) 

<5Q+ 

where C = C(Cfj,,C p ,A,p,5,T). From Lemma r3.ll! we know there 
exists a positive constant C = C{C^ C p , A,p, q, 5, T) for which 

{ue^- c 'rdu\ < I . ) ( 4 (we^- c, ) s dv 



a'Q- 



[a — a' 



for every 0<5<a'<Q;<l and for all < s < q < (p — 1)(2 — p/n). 
Moreover for 5Q~ , since u is a positive superminimizer bounded away 
from zero, we can use Lemma 15.21 to get 



e —^-Q~ ■■ log(ue-?- c ') >\})<C v ,_ . 
Therefore, by Lemma [2. 131 we have 

{ue- p ~ c ') q di^j ^ <C, (6.5) 

where C = C(C M , C p , A,p, q, 8, T). Multiplying with (El gives 
the result 

( f Y /Q 
■f u q dadt <Cessinfw, 

\JSQ- J SQ+ 

for every < q < (p — 1)(2 — p/n), where C = C(C^, C p , A,p, q, S, T). 

□ 



We end this paper by completing the proof of Harnack's inequality 
for parabolic minimizers. This is the first point at which we make use 
of the fact that a minimizer is both a sub- and superminimizer. 

Theorem 6.6. Suppose 1 < p < 00 and assume that the measure fi in 
a geodesic metric space X is doubling with doubling constant C^, and 
the space supports a weak (l,p)-Poincare inequality with constants C p 
and A. Then a parabolic Harnack inequality is valid as follows: Let 
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u > be a parabolic minimizer in Q r C Qt which is locally bounded 
away from zero, and locally bounded. Let < 5 < 1 . Then 

ess sup u < C ess inf u, 
8Q- <5Q + 

where < C < oo and C = C(C M , C p , A,p, 5, T). 

Proof. By assumption it is both a parabolic sub- and super minimizer. 
Hence we may combine Lemma 14.21 with Lemma 16.11 to obtain 

f C \ 1/(p ~ 1} f [ N 
ess sup u < I —r—c — I I -f u p ~ l dv 



< C ess inf u < C ess inf w, 
Where 9 = 6(C^p) and so C = C(C M , C p , A,p, 5, T). □ 
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